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Abstract 

Motivated by recent BABAR measurements of the 7* ^rj and 7* ^rj' transition form 
factors, we estimate two-photon exchange contributions to the corresponding cross sections. By 
using a phenomenological model, based on the vector meson dominance, it is argued that the 
expected contributions are small enough not to effect the BABAR results. As a by product we 
predict Br{r]' fi'^fi") = (1.4 it 0.2) x 10^''. Our results might be useful also in high precision 
calculations of radiative corrections to the Dalitz decays of pseudoscalar mesons. 
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INTRODUCTION 

Recently, the BABAR collaboration reported measurements of the e~^e~ — 777 and 
e'^e^ 7]''-/ cross sections and the corresponding 7* — > 7P transition form factors at a 
center-of-mass energy -y/s = 10.58 GeV. In the asymptotic limit of large photon virtualities 
these form factors are determined by perturbative QCD in terms of the decay constants of 
pseudoscalar mesons. However, the comparison of the BABAR result with the QCD pre- 
dictions is hampered by the fact that the determination of the rj and rj' decay constants 
requires taking into account the mixing between these two states and different phenomeno- 
logical models, cited in give different results. Amusingly, different models are needed 
to reconcile asymptotic QCD predictions for the 1] and rj' transition form factors with the 
BABAR result and neither model can explain the measured ratio of these form factors: 
1.10 ±0.17 versus theoretical predictions in the range from 1.6 to 2.3. This discrepancy, first 
of all, calls for a careful examination of possible sources of unaccounted background. 
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FIG. 1: Two-photon contribution to the e^e — > P7 reaction. 

Two-photon contribution to the e~^e^ P7 reaction, described by the diagrams shown 
in FiglU is one such potential background source. As was mentioned in jl| , this background 
is expected to be very small. But no detailed calculations exist in the literature, to our 
knowledge, to support this conclusion. In this paper, we will try to fill up this gap. 

Naively one can think that this contribution shares the {rrie/Mp) suppression factor, 
inherent to the P e^e" decay, because the graphs in FigH are obtained from e+e~ — *■ 
7*7* — > P by the insertion of a bremsstrahlung photon on the electron line. But this is not 



the case as was realized long ago 



m 



in the context of the P — > e"*"e 7 decay. 



As we see from Figd off-shell P 7*7* amplitude is needed to calculate the two-photon 
contribution to the e~^e^ P7 reaction. Therefore we first discuss this amplitude. 

VECTOR MESON DOMINANCE MODEL 

We define ri{Q) — > 7*(/ci, /i) 7*(A;2, z^) amplitude A^i, through the F^^.^. form factor as 
follows 

A^ij = 2 — , /C2) e^i/crr^i ^2 • (-'-) 

Jit 

The pion decay constant /^r ~ 93 MeV is introduced to make the form factor dimensionless. 
For on-shell_£hotons_and in the chiral limit, = 0, the -Fr?7*7* form factor is fixed by chiral 



anomaly 
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while for large photon virtualities QCD predicts [ld,lll|,ll2,ll3| ~ asymptotic behaviour. 

However, in the interim region the fundamental theories tell little about the form factor 
behaviour and one has to rely on phenomenological models like the Vector Meson Dominance 
(VMD) model 14, Historically this model emerged in attempts to understand photon 
interactions with hadronic matter and in its extreme form asserts that hadronic interactions 
of the photon proceed exclusively through known vector mesons. At low energy region of 
the lightest vector mesons (p, u and 0) VMD was extremely successful in describing a wide 
range of experimental data [16[. This success naturally rises a question whether there is a 
theoretical justification of the VMD from modern perspective of the Standard Model. 

One of the cornerstones of the Standard Model is gauge principle. The naive 7 — p 
direct coupling is not compatible with gauge invariance as it leads to the photon acquiring 
an imaginary mass. Nevertheless it was shown by Kroll, Lee, and Zumino ^3] that the 
complete vector meson dominance is consistent with gauge invariance provided these mesons 
are coupled only to conserved currents. In phenomenological VMD based applications an 
additional terms in the interaction Lagrangian implied by the Kroll-Lee-Zumino analysis 
(for example a photon mass term) are usually neglected as they are of higher order in a. 

Another guiding principle of the Standard Model is symmetry. At low energies we cannot 
solve QCD effectively but its symmetries still guide us in constructing effective theories of col- 
orless hadrons, which are the only relevant QCD degrees of freedom in this non-perturbative 
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region. Pseudoscalar mesons play a special role in this game as they are Goldstone bosons as- 
sociated with the spontaneous breaking of QCD chiral symmetry (in fact would-be Goldstone 
bosons because the chiral symmetry is explicitly broken by nonzero quark masses). Therefore 
their low energy interactions, encoded into a chiral effective Lagrangian, are uniquely de- 
termined from symmetry considerations in terms of a few phenomenological parameters like 
the pion decay constant f-^- The resulting Chiral Perturbation Theory (CHPT) is commonly 
considered nowadays as the effective field theory of the Standard Model at low energies (for 



review see, for example. 
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Although the general method how to add vector mesons (and any other particles) to 
chiral Lagrangians was formulated long ago ^ Is^, usually some dynamical principle is 
needed to reduce number of the free parameters and enhance predictability of the theory. 
Hidden Local Symmetry (HLS) approach [23, 24 1 is considered to be the most convenient 
scheme to deal with vector mesons. Hidden local symmetries were initially discovered in the 
supergravity theories and after realizing that they are a common feature of any nonlinear 
sigma model Bando et al. suggested to consider vector mesons as dynamical gauge bosons 
of hidden local symmetry of nonlinear chiral Lagrangian [2^ There is nothing special 
about this "hidden symmetry". Simply it is a language that makes power counting in 
derivative expansion more convenient when vector mesons are light Q- 

When electromagnetism is introduced in the HLS Lagrangian j26| one finds that generally 
there is a direct coupling of photons to charged pseudoscalars. Only for particular choice of 
parameters, 

( ftrdfrmT \ 1 



(2) 



one recovers the complete vector meson dominance. Equation Q is the celebrated 
Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin (KSRF) relation 27J. Therefore, from this 
perspective, VMD is not fundamental derivative of the Standard Model but rather just a 
lucky dynamical accident of the three flavour QCD 3|. 

Nevertheless the KSRF relation (j21) and hence VMD is experimentally well satisfied. 
Therefore in our estimates we will use VMD form factors. We do not assume SU{3) relations 
between various coupling constants except relative phases, but determine their magnitudes 
from phenomenology. The relevant coupling constants are defined as follows. 

Each photon- vector meson vertex gives a —iegv-yMy factor in the matrix element. The 
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coupling constants gy-y are assumed to be positive and can be determined from the electronic 
widths T{V e+e"): 



viy 

Using the PDG data 



e e 
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, we get 



Myg'v^ 1 + 2 



r Ami Ana^ , ,^ 

1 f ~ ql^My. 



= 0.2014 ±0.0016, = 0.0586 ±0.0010, = 0.0747 ± 0.0012. (3) 
The Tj V^{p)V^{q) transition gives a factor 

This defines dimensionless positive constants g-qvv- Here Sp^u) = 1 and = —1, which 
together with the definition of the gv-y constants corresponds to relative phases expected 
from the SU{3) symmetry with standard mixing angles. Neglecting gr^^i,^, which vanishes 
for the ideal (p — u mixing and nonet symmetry and hence is expected to be small, other 
coupling constants can be determined from the r{V rj'-j) decay widths assuming VMD: 



a 
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The results are 



g^pp = 0.723 ± 0.067, g^^^ = 0.735 ± 0.054, g^^^ = 0.858 ± 0.019. 



(4) 



Now VMD completely determines the F^^.^. form factor in terms of the above given coupling 
constants: 



As an immediate check, one can calculate two photon decay width r(?7 — > 27) using this 
form factor. It is convenient to express the result as the following sum rule 



gvppdp-yMp ^ grjujuj gl^M^ gj^^^g'^^M^ 



+ 



■ (5) 



r(r, - 20O = ^ 



V 



M2 - M2 



(6) 



This relation, which remains valid even for non-zero grjcpuj, is well satisfied experimentally: it 
gives T(r] — > 27) = (0.582 ± 0.085 keV), while the experimental width is 29| T{r] — > 27) = 
(0.510 ± 0.026) keV. 



For r( couplings, we can use r( 
with 



r(^\ V{r( — > 077) and r(?7' p7) as inputs, along 



2 \ 3 



to obtain 



c/^/pp = 0.624 ± 0.044, g^,^^ = 0.724 ± 0.067, ^^/^^ = 0.886 ± 0.059. 



(7) 



In terms of these coupling constants, the F, 



Fri' 



form factor has the form 



+ 



9n'H 9%^4, 



.'rrv'^1,'^2; (fc2_^2)(^2_;\^2) ■ {kl-M^J{kl-Ml) ' {kl - Ml){kl - Mir 
Relative phases are again the ones that follow from the nonet-ansatz SU(3) symmetric 
interaction Lagrangian 

where 

W 

A3, As being the standard Gell-Mann matrices, and the rj — rj' mixing is given by 

r](^s) = cos 9 pT] + sin 9 pT]' , r](^i) = cos 9 pi]' — sin 9 pr], 9p^— 20°, 
while the cj) — ui mixing is assumed to be ideal: 



-^(p°A3 + ^(8)A8) + ^a;(i), P= -^r/(8)A8 + -^r/(i), 



^(8) = v^'^ + v^^' 
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Sum rule that follows 

r(V ^ 27) 



3M3 
2a 



Mr,' 



M^,.-Ml\ 



r(r/' ^ P7)r(p 



+ 



(9) 



r(r7' — > co'7) r(ct; 



+ 



3r(0->V7)r(0 



M2 - 



is well satisfied and supports our assumptions. It gives r(?7' 27) = (4.56 ± 0.32) keV and 
the experimental width is |29| 1(77' -> 27) = (4.29 ±0.15) keV. 

But r(P ^ 27) widths check the Fp^*^* form factors only for on-shell photons and 
the two photon contributions into e+e^ P7 we are interested in depend on the off-shell 
behaviour of these form factors. Note that © and (jH]) apparently violate the QCD prescribed 
~ asymptotic behaviour when both photon virtualities are large. Therefore now we 
turn to the r] /i+p~ decay, where the off-shell behaviour of the form factor does matter, 
to demonstrate that (0) still gives reasonable estimate. 



fi+fi- DECAY 



There is an intensive literature devoted to the rare decays of 



_2seudoscalar mesons into 



a lepton pair (for a review and references see, for example, |3C . l3ll . l32l . l33L l34|). Detailed 
pedagogical calculation in the framework of VMD is given in 35|. In somewhat different 
manner, this calculation is reproduced in js^. We follow these references in spirit but differ 
in technical details. 

The leading contribution to the ri{Q) fi'^iQ — p, s+) + fi~{p, s_) decay comes from the 
diagram shown in FigO 




FIG. 2: The leading contribution to the r] fi^fj, decay. 



The corresponding Feynman amplitude is 



4 4 

6 e 
^ = — u{p, sJ)Mv{Q - p, s+) = 



where 



M 



v{Q - p, s+)u{p, s^)M 
dk -f^{p -k + m/,)7, e'^'^'^k^Qr „ 



(27r)4 [{p - ky - ml] fc2 (Q - k) 



2 '77 7' 



{k\{Q-kY). 



(10) 



The lepton pair from this decay has total angular momentum J = and hence is either 
in singlet ^5*0 or in triplet ^Pq state. But the triplet state has CP-parity (—1)'*+^ = 1 
which does not matches the negative CP-parity of the t] meson. Therefore, assuming CP 
invariance, v{Q — p, s+)u{p, s_) in (fTUI) can be replaced by the projection operator to the 
singlet state for the outgoing system 

1 _ 
V{Q -p,p) = -j=[v{Q - p, +)u{p, -) + v{Q - p, -)u{p, +)] . 



This projection operator was calculated in j3l| with the result 



V{p+,p- 



2V2t 



-2m^ {p+ + p_)75 + - e^,^r{P-Pl - plP-)cr^'' + t^^ 
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where t = (p_ + and a^'^ = ^{j^'j'^ — 7^7^)- Therefore 

^ dk e'''"'^KQrL 



-Pr;7*7* {k , {Q k) 



with 



But then 



L^,u = Sp V'j^ip -k + m^)7^ = -2iV2 ^ e^y^^k^Q' 



In the above expressions REDUCE was used to perform Dirac algebra and calculate 
traces. 

Finally, we get in the standard way 

2 



2a , , / m,, 



TT 



-M 



4"^?, , ,9 

1 - ^ LR 



(12) 



where j38, 



39| 



7r2 7 M2 P (Q _ [(p _ _ ^2 

Remembering VMD expression (jSJ for the F^j^*^* form factor, we can write 



(13) 



R = gr^pp I{Ml M^) + g^^^ gl^ I{Ml Ml) - g^^^ gl^ I{Ml Mf 



where the master integral has the form 



(14) 



ml Ml 



dk 



M^k^-{k-Qy 



J M2 (g - ky [{p - ky - mi] (p - md [(q - - m|] 



(15) 
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The imaginary part of this integral can be calculated by using the Cutkosky rules 
4l|. When intermediate masses Mi^2 are greater than the //-meson mass, only two-photon 
cut contributes to the discontinuity of /(M2,M|) (which is twice its imaginary part) and, 
therefore. 



disc I {M^, Ml 
where 



-2m] 



dk{~k^) 



5^{k^)54{Q - kf 



Ml Ml 



k^-2p-k {k^ - Mf) [{Q - ky - Mi 



21 ' 



6+{e) = Q{ko)6{e) = —S{ko-\k\) 
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Because of (5-functions, calculation is straightforward and gives the well known model- 
independent result 



ImI{M,,M,) = ]^discI{MlMl) = ^ ln[^, P = ^ 



4m?, 

1 ^. 

M2 

rj 



(16) 



Calculation of the real part is much more tricky and as a first step involves the following 
algebraic identity 



[k'^Q^ - (k ■ QY] Ml Ml _ \{Mf, Ml 



(g2 _ M|)2 _ r 1 1 



D,D 



2^3 



1 1 



(17) 



where 



Di = D2 = {Q- kf, D3 = P-2p- k, Di = e- Ml, D5 = {Q - kf - M; 



2 ' 



and A(x, y, z) = + y'^ + — 2xy — 2xz — 2yz is the triangle function. Using this identity 
and introducing the dimensionless variables 



mr 



M| 



we get 



ri 



/(M^M|) = J(n,r2) + J(0,0)-J(n,0)-J(0,r2)--^[/i(p2)-/i(0)]-5[/2(pi)-/2(0)] 



(18) 



with 



/l(P2) 



7r2 y [A;2 - 2p ■ A;] [(Q - A;)2 - M| 



21 ' 



/2(Pl) 



dk 



21 ' 



and 



1 ^ /" 
^(ri,r2) = -A(l,ri,r2)fif(ri,r2), fi'(ri,r2) = r / ■ 

4 TT^ J 



[A;2 - 2p ■ A;] [A;2 - Mf] [(Q - A;) - ...2 
Integrals with two denominators (/i and /2) are easy to calculate by using 

1 



1 



dx 



AB J [xA+{l-x)B]^ 

and the dimensionally regularized integral {'Je is the Euler constant) 
I r dk r(e/2) /2 



71^ J [A;2-A]2 (7rA)^/2 



— 7^; — Invr — /nA , e = 4 — (i— >0 
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The result is l35 



36| 



/i(p2)-/i(0)= dxln 



and 



m^x^ + MKl-x) 1 



2p2 



/2(pi)-/2(0) 



2pi 



Inpi + Jl - 4pi In 



lnp2 + \/l - 4p2 In 



1 + 



1 + yi - 4p2 
1 - VI - 4p2 



(19) 



1-^1^4;^ 



(20) 



As for the integral with three denominators g{ri,r2), we use Feynman parameterization 



1 



1 X 



2 / dx 



dy 







[(x-y)Ds + {l-x)D^ + yD,]^ 



and the integral 



to get 



1 X 



71^ J [P-A]3 2A 

dy 





Let us shift the y- variable as follows: 



r (x^ + y2) — + ?/ + rix — r2y — ri 



Then 



with 



But 



y y + ax, a = 

1 {l—a)x 

g{ri,r2) = dx 



l + P' 

dy 



—ax 



f{y)x + h{y) 



fiy) =ri + ^-p^ (r2 - 1) - /3y, h{y) = -ry"^ + (1 - r2)y - ri. 



1 {l-a)x 1 {l-a)x 1 -ax 1-a 1 



-a 1 



y (ix y dy = j dx j dy — j dx j dy = j dy j dx — j dy j dx 

-ax y/(l-a) -y/a 

and after performing the trivial x-integration we end up with 

2/{l+/3) 



C/(n,^2) 



dy 

W) 



(l-/3)/(l+/3) 



/ 



dy 

m 



In [/(y) + %)]-ln 



In [/(y) + %)]-ln 



f{y)y{l + f3) 



f{y)yil + P) 



h{y) 



+ h{y) 



(21) 
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We need only the real part and, therefore, integrals in (|7T|l are of the type (a, a, b, A, B, C 
are some real constants here) 

rv 

dy 



F{a; a, b; A, B, C) = Re f — ^ In {Ay^ + By + C). 

J ay + b 



(22) 



A particular case of such type of integral 

dy 

In {y - yi){y - 1/2) - In {yo - yi){yo - ^2) 

y-yo[ 

was considered in detail in (see also Using their result 







ReZ = Re 



yo 



Li') 



yo - 1 



yo 



Lin 



\yo-yij \yo-yij \y0-y2, 

we can readily calculate ()22p in terms of dilogarithms: 



yo - 1 
.yo - y2, 



F(a;a,b;A,B,C) = -Re \ \n(l + ay) In (a^-B- + c] + 
a [ \ b J \ a'' a J 



Lin 



Li2 I ^_lL^ I _^ 



where 



b + yia ) \b + yia 

-B + - AAC 



b + y2a 



— Lin 



b + aa 
b + y2a^ 



(23) 



yi 



2A 



1/2 



-B - - 4AC 



2A 



Returning to ()2H1 . we get finally 

2 



Reg{ri,r2) 



F 
F 
F 
F 



,l + /3' 
,l + /3' 
1 + /3' 





l + P 


(?^2 


-1); 


-r, 1 - /3 - 


P,n + 


1-/3 
l + P 


(^2 


-1); 


-r, 1 - /3 - 


P,n + 


1-/3 
1 + P 


ir2 


-1); 


(1 + /?)' 
4 




1-/3 
1 + /3 


(^2 


-1); 


4(1-/3) 



./3 \ 
— (1 + ri -r2),-ri + 



■ri, -ri 



(24) 



Some arguments of dilogarithms in (j24j) are complex. Therefore for numerical evaluation we 
need an algorithm to evaluate this function for complex argument. We use the algorithm 
described in ^J] (see Q] for another algorithm). 

First of all the argument of the dilogarithm is brought in the region |x| < 1, —1 < 
Re{x) < I by using the functional identities 



Li2(x) 



IT 



Li2(l - x) + — In (x) In (1 - x), Li2(x) 
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Then the series expansion 



n=0 



{n+l)V 

truncated at n = 20, is used with Bn Bernoulh numbers and z = — In (1 — x). 

Now we have all ingredients at hand to calculate numerically r(r7 — > fi^fx^). The result 

is 

Br{7] /iV) = ^^^"1^^^ ^ = (5.2 ± 1.2) X 10"^ (25) 



The quoted uncertainty is dominated by uncertainty in the Qnpp coefficient. The experimental 
number is Br{ri — > yU+yU") = (5.8 ± 0.8) x 10^^ and again we observe a good agreement. 
However, it should be mentioned that the dominant contribution into (f^3j) is given by the 
imaginary part of the amplitude which is in fact model independent: 

Re{R) ^ -0.015, Im{R) ^ -0.074. 

To check our formulas and computer code, we have make sure that they reproduce numerical 
values of 2/(M^, M^) and of the analogous integral for the vr" — > e~^e~ decay given in js^. 

Analogous considerations apply to the rj' decay if we use the form factor dHJ 

with coupling constants But now the imaginary part of the amplitude has an additional 
contribution from the on-shell pj and intermediate states. This contribution is calculated 
by using the Cutkosky rules with the result 



TT 



AImIiMt,Mi) = - — 



1 - 



Ml 
Ml 



1 - 



Ml 

Ml 



In 



l + P' 



(26) 



The numerical calculation yields the branching ratio 



Br{r]' fi+p-) = (1.4 ± 0.2) x 10 



-7 



(27) 



Now the imaginary part of the amplitude is only 1.5-times larger in magnitude than the real 
part: 

Re{R) ^ 0.064, Im{R) ^ -0.093. 

Unfortunately there is no experimental number for the rj' ii~ branching ratio to com- 
pare with (jTfj) . 
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TWO-PHOTON CONTRIBUTION TO e+e" ^ r/7: CONSTANT FORM-FACTOR 



The amplitude of the e^{p+, s+) + e s_) r7(g) + 7(fc, e) process is given by 

A = v{p+, s+)T''{p+,p_, k)u{p-, S-) e*, 

where gauge symmetry and invariance under parity and charge conjugation dictate the 
following decomposition of in terms of four independent invariant form factors (up to 
irrelevant additional structures proportional to k^) 45[ 

= P(x-, X+) [{k ■ p+) p'l-{k-p_) p1] 75 - ^ T(x_, X+) ^^''k,j5+ 



A+{x~,X+) - ■P+)7'' 75 - kp--{k-p^)-i 



75, 



with 



(28) 



P(X-,X+)=nX+,X-), ^±(X-,X+) = ^t(x+,X-), T{x-,X+) = T{x+,X-)- (29) 

Here x~ = '^k ■ p_ and x+ = 2 A; ■ p+. 

The invariant form factors P, A± and T can be projected out from by a suitable 
projector operator A^, F = P, A±, T according to the the formula 



F = Sp [A^ {p+ - nie) (p_ + me)] 



(30) 



We are interested in the mg — >■ limit and in this 
expressions differ somewhat from ones given in 

1 



imit these projectors have the form (our 



where 



2A252 



-i e'""'^k^p+„p^r 



2 A'^75 



^± 16(A2)2 



6'^ k ■ p^ k ■ p- 
2A2 ± p"^ k ■ {p t 5) k ■ 5] \k ■ pY - kp'' 



± 



p' (2^ - k ■ {pT S) k ■ p^ {k ■ 6^ - kd^)"^ 



A^ 



4A- 



(31) 



^ = {k ■ P+)p- - {k ■ pS)p+, 5 = p^-p^, p = p++p^. 
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, it can be 



Using these projector operators and calculating traces by means of REDUCE 
found that P and T form factors vanish in the mg — > limit for diagrams shown in FigO 
Therefore, in this limit we get in the standard way 



(e+e- r/7) _ M^^ 



l + cos0)2|A+|2 + (l-cose)2|A^|^ , (32) 



dn A^TT^ \ s J 

where a^'^"'^ (e^e~ — >■ r/7) denotes purely two-photon contribution into the e+e~ — >■ 777 cross 
section, s = (p_ + and 9 is the flight angle of the photon with respect to the electron 
momentum p_ in the center of mass frame. 

To calculate (jH^ . we need to find out only the y4_(x_,x+) form factor, because of the 
charge conjugation symmetry relations ()29|) . The following method [3j for the form factor 
determination seems to be more convenient than the direct use of the suitable projector. 

In the massless electron limit 

r'^ = [Ai p>t - A2 pI] A; 75 + ^3 ^75, (33) 
where Ai = —A_, A2 = —A+ and A3 = k ■ A2 — k ■ A2. Using the identity 

it can be shown that the first two diagrams from Fig^ contribute to A^ alone. Therefore 
to determine A_ = —Ai it is sufficient to consider only the third diagram with photon 
emission from the internal line. The contribution of this diagram in the massless electron 
limit looks like 

r(c) f 7.(P+ - 07/.(P+ - i - khx e'^'^n^q, „ / _ ^^2^ .o.^ 

/.i (27r)4 /2 _ /)2 _ /)2 _ / _ ^.7*7*1', W UJ- K^^) 

As follows from the coefficient of the k^p_^ term in Sp{'y^p_T^'y^p^) is — 4p_ -p+Ai = 
'^P- ■ P+ A_. Besides, (jMj) is free from ultraviolet divergences even for constant form factor 
Frfy*'y* ( polnt-likc T]). Thcrcforc in this section we assume, as in [3], that the -F^-y*^* is just 
a constant. Then, to find out the A^ form factor, we combine Di = D2 = {q — lY, 
-D3 = {p+ — 0^ ^^'i -^4 — {P+ — ^ — ky denominators in ()H4|1 by using 

1 } } 1 

= 3! / dx-i ■ ■ ■ dxd 5 fy^ Xj — 1) 7-, 

DiD2D3D^ i i ^ \^ ' ) [xiDi+X2D2 + x^D^ + XiDiY' 

shift the variables in the momentum integral according to 

1^1 + X2P- + (1 - Xi)p^ - {X2 + X4)k, 

14 



calculate the trace S'p(7y]3_r(f)75p+) with the help of REDUCE [37|, replace in resulting 
momentum integrals 



and separate the coefficient of the k^p^^^ term. As a result we obtain 

1 



J dxi ■ ■ ■ j dxi 5 Xi — 1 



X 



(2vr) 



dl 





[1 — 4X2 — 2X4) — 25X2X3 + 2X+(X2 + X4 — 1)(X2X3 — X1X4) 
[P - X2X3X- - X1X4X+ + ^« X1X2] - 



(35) 



Our result (jH^j) differs somewhat from the corresponding expression in j3[ . For this reason we 
have cross-checked it by an independent calculation where part of the above given program 
was performed manually. 
Momentum integrals read 

I i f „ 1 1 



6A^' 



Therefore 



± 1 
J dxi ■ ■ ■ j dxi 5 Xi — 1 



X 



1 — 4X2 — 2X4 _|_ SX2X3 — X+(X2 + X4 — 1)(X2X3 — X1X4) 



A 



A^ 



where 



A = X2X3X- + X1X4X+ - M^xiX2. 
Let us introduce dimensionless variables 



and note that 



Then 



with 



ri Tj 



19 1 X1X4 19 1 X2X3 



2 + (X+ + X_ - 1; 



d 



/2(X+,X_) + 



(36) 



(37) 
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2 + (l-X, 



d 



d 



and 



2.(2) _ 



[/2(X+,X_)+/4(X+,X_)] 



dX, 



Here the master integrals are 

1 1 

/o(x+, x_) = y (ixi ■ ■ ■ y ^^2:4 5 Xj — 1 



1 1 

In{X+,X_)= dxi--- dxi5{Vxi-l. 

J J ^ X1X2 + 1 — A_ X2X3 + 1 — A+ X1X4 

The master integrals (|HS|) can be considered as integrals over a tetrahedron shaped 3- 
dimensional domain in the (xi, X2, X3)-space. By means of variable transformation [3| 



X1X2 + (1 — X_)X2X3 + (1 — Xj^)xiXi 



Xr, 



(38) 



xi = xz, X2 = y{l-z), X'i = z{l-x), X4 = (1 -z)(l -y), 



(39) 



with the Jacobian 



d{xi,X2,X'i) 



z(l-z) 



d{x,y,z) 

the integration domain transforms into a unit cube in the (x, y, z)-space and we obtain 



111 
h = J dx J dy J dz 
000 

1 



y(i-^) 



;i - x+)x + (1 - x_)y + (x_ + x+ - i)xy 



1 f In [1 - X+ + X+i/ + ie] - In [(1 - X_)y + ie] 

2/ ^ i-x++(x_+x+-i)y y- 



(40) 



Analogously 



1 _1 f ln[l -X_ +X_x + ze] -ln[(l -X+)x + ze] 
^^ + ^^-2^o-^Jdx i-X_ + (X_+X+-l)x 



(41) 



In the above formulas the ie prescription for propagators was explicitly restored. 

First of all let us calculate the imaginary part of the amplitude which is in fact model 
independent as far as it is dominated by the on-shell two-photon intermediate state. In ()4U|) 
and ()41|) the imaginary parts originate from negative arguments of the logarithms. Therefore 

1 

Im {I2 + h)-- 



TC 

2 



dx 



1-1/X- 



1 - X_ + (X„ + X+ - l)x 
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TT 



In- 



and 



Im lo 



2(X_+X+-1) (X_-1)(X+-1)^ 

ydy 

1 - X+ + (X. + X+ - l)y 



TT 

2 



TT 



2X, 



1-1/X+ 



+ 



X+(X+ - 1) 



In 



X_X4 



X+ + X_-l (X+ + X_-l)2 (X_-1)(X+-1 



Inserting these expressions into (jHTjl we get finally 



Im A_ 



777*7* 



X, -x_ 



8nfM X_X+(X_-1) 



cos^ 



1 — COS0 



MS 



A/2 



and the corresponding contribution into a^'^'^^ (e"^e — > r/7): 



(27) 



to 



3 /M^V 



2\ 3 



M2 

77 



r(7^ ^ 77) 

Mr, 



cos^ 



cos^ 



(42) 



2 ' 



(43) 



where the constant F^^y*^* was expressed trough the two-photon width T{ri ^ 77). Integrat- 
ing in the limits 30° < 9 < 150°, and assuming s = 112 GeV^ we find numerically 

af^^(e+e-^r/7)^3.4-10-^fb, ^^(e+e" ^7) ^ 4.5 ■ 10"^ fb. (44) 

These numbers are quite small compared to the reported BABAR cross sections Q| 

(T(e+e- ^ 777) = 4.5l{:? ± 0.3 fb, a(e+e" ^ ri'-f) = 5.4 ± 0.8 ± 0.3 fb. 

But the amplitude has also a real part which in fact dominates for the constant F^^*^. form 
factor. 

The real part of the integral PT|) can be evaluated by using a substitution t = 1 — X_ -|- 
(X_ + X+ — 1) X and the integral 

B . 



In (1 - Ct) 



t 



dt = Li2{CA) - Li2{CB), 



A,B and C being some constants. As a result, we get (note that Li2{l) = tt^/G) 



Reilo + h) = Re 



2 



2(X_ + X+ - 1) 
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By using the Abel identity 
In (1 — x) In {I — y) = Li2 



46 



471 ] (still valid for real parts if x > 1 and y > 1] 



X 



1 — X 



■ Li2{x) - Li2{y) - Li2 



xy 



the above expression can be rewritten as 



Re {h + h) 



ReG{X_,X+] 
2(X_ 



:i-x)il-y)J' 
(45) 



where 



TT 



G{X_,X+) = In (1 - X+) In (1 - X_) + Lt2{X+) + Lt2{X_) - —. 

D 

The real part of the integral (jin|) can be calculated analogously with the result 

X+ - 1 In (X+ - 1) (X+ - 1) G{X^, X^ 



Reh 



2(X„+X+-1) 
Using (gSI) and (jUj), we find 



Re 



In 



X_ - 1 



X4 



(46) 



(47) 



Re A 



(1) 



G(X_,X+) 



4fc ■ p_ 



in agreement with . But according to our results A^"* is only a part of the A_ form factor, 
for which we obtain 



ReA_ 



e^K 



X+ - 1 In (X_ - 1) In (X+ - 1) 1 



In contrast to [3|, all dilogarithms have canceled out. 

Having at hand ReA^, we can find the corresponding contribution in the cx^^'^-' (e"*"e 
777) cross section: 



(4J 



(icr 



(27) 



■q 



47r2M2 



|F(cos^)|2 + |F(-cos( 



(49) 



where 



and 



F(cos( 



In 



1 + cos^ ln(X_-l) ln(X+-l^ 



cos 9 



X4 



1 + 



X_ 



X4 



1 
2 



2M2 



'1 ± cos I 



Numerically, for 30° < ^ < 150° and s = 112 GeV^ we find 



a 



(27), 

R ' 



e+e- r]-f) ^ 0.5 fb, ap\e+e- t/'t) ~ 0.8 fb. 



(50) 



(51) 



These cross sections are quite substantial and they could cause problems with the BABAR 
analysis if the constant F^^^*^* form factor would not be unphysical. In reality, however, 
the form factor drops quickly at large photon virtualities. Therefore in the next section we 
consider more realistic VMD form factor. 
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TWO-PHOTON CONTRIBUTION TO e+e" ^ 7/7: VECTOR MESON DOMI- 
NANCE MODEL 

The VMD -F^-y*^* form factor is given by (0). Its insertion into (j34j) produces momentum 
integrals with six denominators. Their treatment becomes less formidable if we use the 
identity 

1111 



DiD2D:,D^D^D^ D^D^D^D^ D^D^D^D^ D^D-^D^De D^D^D^D^ 



with 



Then the procedure described in the previous section gives 



with 



Here 



i{Ml Ml) = J{Ml Ml) + J(0, 0) - J{Ml 0) - J(0, Ml). 

d 



-lf~M\ 



(52) 



(53) 



M^J{Mt,M^) 



2 + (1 - X, 



d 



2 + {X+ + X_- 1) 

d ^ 



/2(X+,X_ 



1 + (1-X 



dX^ dX+ 
' d 



dX_ dX, 



dX_ 

/2(X+,X_)+/4(X+,X_ 

^ ' /o(X+,X_), 



(54) 



and the new master integrals are 

/o(X+,X_) 



1 1 
= j dxi ■ ■ ■ y 



1 1 
/n(X+,X_) = / dxi ■ ■ ■ dx4 



HEx^-l) 



X1X2 + (1 - XJ)X2X2, + (1 - X+)xiX4^ - XiTi - X2r2 
SiX^i- 1) Xn 



■ (55) 



X1X2 + (1 - X_)X2X3 + (1 - X+)xiX4 - XiTi - X2r2 
u u 

The presence of the additional Xiri +0:2^2 term in the denominators compared to ()38|1 makes 
the evaluation of these integrals, of course, more complicated but not substantially different 
from the evaluation of (j38|) . Therefore we present only final results. 



Re {12 + h) 



(X_+X+-l) 



Re I dz{l- z)G{X'__,X' 



19 



Relo 



Relo 



(x_ + x+-i)- 

1 



Re / dz{l-z) 



In 



x; - 1 In (x; - 1) (x; - 1) G{x'_,x'_^ 



XL - 1 



X\ 



+ 



X'_ + XL-1 



Re / dzG{X'_,X' ). 



(X_+X+-l) 
Here the function G{x, y) is defined by ()46p and 

2 [ri - (1 - z)X_ 



(56) 



X' 



xi 



(1 - z) [r2 - zX+] 



(57) 



ri2; + r2 (1 — -z) — 2;(1 — z)' ^ riz + r2 (1 — 2) — z(l — z) 

We need only the real parts because the imaginary part of the amplitude is dominated by 
the two-photon cut and, therefore, is the same as for the constant Fri^y*^* form factor. 
Using (jSEI), we get after some algebra 

M^Re J{Ml, Ml) = ^Re J dz{l - 2z) G{X'_,X'^) + 



Re / dz 



z{l - z) 



r2 + z{X_ - 1) 



ri+X+-l 



ri + (X+-l)(l-z; 



L(xi,x;)-i 



where 



Note that 



and therefore 



Lix,y) = In^ + ^"^^-^^ - h^y^^ 
X — 1 X y 



Ml h + z(X„ - 1)] 



2M| / 



+ z\l ^ (l-cos( 



ReJ{MlMl) = -Rej{n,r2), 
s 

where the dimensionless function j(ri, is given by 
1 } 

ra) = -j dz{l - 2z) G{XL, X;) + 



(58) 



dz 



z{l - z) 



+ z(l-^) (1-cos^) 



ri + X+ - 1 



ri + (X+-l)(l-^) 



L(x:,x;)-i 



The corresponding contribution in the o"*-^''''' (e^e — > r77) cross section looks like 



d'y^R^\e'^e — > 777) 



M 

X _ 



2\ 3 



|/(cos^)|2 + |/(-cos( 



(59) 



(60) 
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where 

/(cos^^) = (1 -cos^) [gr,ppg%i{MlMl) + g^^^gl^i{MlMl) - g^^^gl^i{MlMl) 
and 

i{Ml Ml) = Jin, r^) + j(0, 0) - Jin, 0) - j(0, r2). 

Now we have all ingredients to find the cross section under interest numerically. The results 
are (for the same kinematic conditions as before) 



P\e+e- ^r]-f) ^1.6-10-^ fb, ag^^(e+e- ^ r/'7) ^ 1.2 ■ 10"^ fb. (61) 



For e~^e~ rj'y the real part of the amplitude still dominates although not as drastically as 
for the constant F^^*^* form factor. For e^e" —>■ Vj^ the real part is even somewhat smaller 
than the imaginary part. 

CONCLUSIONS 

In this work we used a phenomenological, VMD inspired F^^.^* form factor to estimate 
the expected contributions into the e^e~ —>■ rij and e~^e~ —>■ amplitudes from the two- 
photon exchange diagrams of Fig^ The results are given in (j44|) and (jblj) . The numbers 
obtained are too small to be of any importance for the recent BABAR analysis but 
indicate that interference effects of the order of several percents are expected in future high 
precision studies of corresponding transition form factors. 

Our results might be of some relevance for accurate evaluation of radiative corrections 



to Dalitz decay P ^ e^e 7 



especially in light of some discrepancy with earlier studies 



but we do not pursue this line of thought further in this article. 
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